Interedge magnetic coupling in transition-metal terminated graphene nanoribbons 
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The magnetic structures and interedge magnetic couplings of Fe, Co and Ni transition-metal terminated 
graphene nanoribbons with zigzag (ZGNR) and armchair (AGNR) edges are studied by first-principles cal- 
culations. Fe-ZGNR is found to show antiferromagnetic (AF) coupling between two edges, while the interedge 
coupling of Co-ZGNR is ferromagnetic (FM). For Fe-AGNRs and Co-AGNRs, increasing the interedge distance 
we follow oscillatory transitions from FM to AF coupling with a period of about 3.7 A. The damped oscillatory 
behavior indicates a Ruderman-Kittel-Kasuya-Yosida type interedge magnetic coupling and the oscillation pe- 
riod is determined by the critical spanning vector which connects two inequivalent Dirac points in the graphene 
Brillouin zone. The two edges in Ni-ZGNR are decoupled independent of the ribbon width and Ni-AGNRs are 
found to be nonmangetic. 

PACS numbers: 73.22.Pr, 75.30.Et, 75.75.-c 



Graphene, a monolayer of carbon atoms packed into a hon- 
eycomb lattice, continues to attract immense interest, mostly 
because of its two-dimensional stability, unique band struc- 
ture and other unusual physical properties 1 . In particular, 
cutting graphene along two high-symmetry crystallographic 
directions produces quasi-one-dimensional periodic strips of 
graphene with armchair or zigzag edges, usually referred to 
as graphene nanoribbons (GNRs). The zigzag edge GNR 
(ZGNR) is theoretically predicted to be magnetic with two 
spin-polarized edge states, which are ferromagnetically or- 
dered but antiferromagnetically coupled to each other through 
the graphene backbone^, while the armchair edge GNR 
(AGNR) is found to be nonmagnetic. The interedge magnetic 
coupling in ZGNRs has attracted considerable attention 3-6 . 
Antiferromagnetic coupling of the two zigzag edges in ZGNR 
can be explained in terms of interactions between the mag- 
netic tails of the edge states 3 , since the C atoms always stand 
at the opposite sublattices of GNR at the two zigzag edges. 
The magnitude of the interedge magnetic coupling shows a 
w~ 2 dependence as a function of the ribbon width w^&. 

Most of the studies of interedge magnetic coupling on 
GNRs focus on ribbons with zigzag edges and hydrogen ter- 
minations. We have previously reported that metal terminated 
GNRs can also exhibit magnetics behavior 7 . Furthermore, 
GNRs with ferromagnetically coupled edges terminated with 
transition metals show high degree of spin polarization at the 
Fermi energy, and thus can be excellent candidate for spin- 
tronic applications. The interedge magnetic coupling in Fe 
terminated ZGNRs has also been studied by Ong et ai— very 
recently, showing that the coupling is antiferromagnetic and 
the strength decreases with increasing ribbon width. 

For a full understanding of the interedge magnetic cou- 
pling in TM-GNR systems, in this paper, we present a first- 
principles study of zigzag and armchair GNRs terminated 
with Fe, Co and Ni 3d transition metals, focusing on the width 
dependence of the magnetic coupling between two edges. In- 
terestingly, the behavior of interedge magnetic coupling is 
found to differ significantly with different type of metal ter- 
minations or ribbon edges. We also find a damped oscilla- 
tory behavior of interedge magnetic coupling and the oscilla- 
tion period is determined by the critical spanning vector which 
connects two Dirac points in the graphene Brillouin zone. 
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FIG. 1: (Color online). The top view of structures of TM-ZGNR 
with N z = 8 (a) and TM-AGNR with N a = 14 (b). The solid (red) 
dots represent transition-metal atoms passivating the ribbon edges. 
The dashed rectangle denotes the unit cell in calculation, (c) and 
(d) are corresponding side views of the TM-ZGNR and TM-AGNR, 
respectively, a, b and L are the two metal-metal bond lengths at each 
edge and the the length of unit cell, respectively. 



The electronic structure calculations are performed using 
density functional theory implemented in the plane-wave- 
basis-set Vienna ab initio simulation package (VASP) 9 . Each 
ribbon is simulated within a supercell geometry containing 4 
metal atoms at two edges, as shown in Fig. Q] A large vac- 
uum spacing of 15 A is used between two edges and between 
two graphene planes to prevent interaction between adjacent 
images. Projector augmented wave (PAW) potentials with 
kinetic energy cutoff of 500 eV are employed in all simula- 
tions. For the exchange and correlation functional we use the 
Perdew-Burke-Ernzerh generalized gradient approximation 10 . 
Brillouin-Zone sampling is done on a grid of 36 x 1 x 1 
Monkhorst-Pack 1 1 ^-points along the periodic direction of the 
ribbon for ZGNRs and 40 x 1 x 1 for AGNRs. The Gaussian 
smearing method is used to treat partial occupancies, and the 
width of smearing is chosen to be 0. 1 eV for geometry relax- 
ations. The geometries are optimized until all forces on all 
ions fall below the threshold value of 0.01 eV/A. To obtain 
accurate magnetic configurations and ensure high accuracy in 
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TABLE I: Structural, energetic and magnetic properties of metal-terminated GNRs. Columns show: metal termination, metal-metal (a or b) 
and metal-carbon (du-c) bond length (in A), the binding energy Et and the formation energy of the metal chain fifijjJLj.. (in eV per metal 
atom), total energy difference between AFE and AF (in meV per unit cell) states, total energy difference between FM and AF (in meV per unit 
cell) states, total magnetic moment m tot of the ribbon in the FM state (in Bohr magneton ^ per edge termination), local magnetic moment of 
the metal atom /i M and its nearest-neighbor C atom yu c (in 
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the calculated total energies of relaxed structures, calculations 
are conducted using Gaussian smearing of 0.05 eV and PAW 
pseudopotentials for Fe and Co which treat the 3p semi-core 
states as valence states. 

Two typical structures of metal-terminated GNRs with 8 
carbon zigzag chains and 14 carbon dimer lines across the rib- 
bon width are shown in Fig.[T](a) and (b), respectively. Here- 
after we refer to a metal-terminated GNR with N a dimer lines 
as a AVTM-AGNR and that with N z zigzag chains as a N z - 
TM-ZGNR where TM stands for Fe, Co or Ni. Both edges 
of the ribbon have the same configuration for each case. We 
find the lowest-energy structures of the ribbon by structural 
relaxation calculations for three possible termination config- 
urations: the linear type (a = b — L/2), the dimerized linear 
type (a + b and a + b = L), and the zigzag type (a + b > L). 
The calculated ground state structural properties are listed in 
Table U for 10-TM-ZGNRs and 17-TM-AGNRs. The results 
show that the most favorable structures for all studied TM- 
ZGNRs are all linear while the TM-AGNRs are all zigzag type 
mainly because of a relatively smaller length L of the AGNR 
unit cell. The Fe and Co terminations dimerize at the edge of 
ribbon due to the Peierls distortion 12 . Changing the width of 
the ribbon will have negligible effect on the structure of the 
edges. The binding energies Eb of the metal atom, defined as 
E h = (£ T m-zgnr - £zgnr)/4 - E^ m , and the formation en- 
ergy of the metal chain, £f h ™ ation = £^°/ 2 - E tm> ^ ^ so 
listed in Table|T] The difference between Eb and E^^^. rep- 
resents a direct binding between the TM and carbon atoms at 
the ribbon edge. It is clearly shown that for all cases the metal 
atoms bond strongly with edge carbon atoms. 

Next we examine the magnetic structures and couplings 
between the magnetic moments in the ribbon edges. Three 
states with different spin configurations of the TM termina- 
tions are considered: (i) antiferromagnetically ordered spins 
at each edge of the ribbon, denoted by AFE, (ii) ferromag- 
netically ordered spins along both edges with the same spin 
direction, denoted by FM, and (iii) ferromagnetically ordered 
spins at each edge with the opposite spin directions between 
the edges, denoted by AF. Total energy calculations are per- 
formed to decide the ground states of the magnetic structures. 

Except for Ni-AGNR which is found to be nonmagnetic, 
each ribbon considered shows spin polarized edges with fer- 
romagnetic ordering of the metal atoms at each edge for the 



lowest-energy state. In Table[TJwe show examples of 10-TM- 
ZGNRs and 17-TM-AGNRs. Besides that the AFE state of 
Ni-ZGNR is found to be not stable at all, it is clearly shown 
that AFE state with antiferromagnetically ordered spins at 
each edge are not favored, with a large energy difference as 
compared to FM or AF states. Though in the Fe-AGNRs and 
Co-AGNRs the metal terminations alternatively bond to the 
opposite sublattices of GNR along each armchair edge, the 
favorable FM (or AF) state is consistent with the ferromag- 
netic ordering found in the ground states for Fe, Co and Ni 
monatomic chains^. In Fe-ZGNRs and Co-ZGNRs the en- 
ergetic disadvantage of AFE state is very obvious, and can 
be explained by the interactions between the magnetic tails 
of the edge spins 3 , since the metal atoms always bond to the 
same sublattices of GNR at each zigzag edges. Generally the 
magnetic moment of the ribbon comes mostly from the metal 
atoms and their nearest-neighbor C atoms at the edges, and the 
edge C atom presents magnetization antiparallel to the nearby 
metal atom, as shown in Table Q] for ribbons with FM states. 
In the AF state the moments at two edges have exact the same 
values but with opposite signs, thus the net magnetic moment 
of the ribbon is zero. 

The magnetic order between the two edges for the ground 
state can be either FM or AF favored depending on the in- 
teredge magnetic coupling. Flipping the spin moments of one 
of the edges will result a total energy change in the system 
for comparable size of the magnetic moments at the metal ter- 
minations. The interedge magnetic interaction strength can 
be identified by the total energy difference A£p^ AF between 
the FM and AF states. Our calculations show that the ground 
states of Fe-ZGNRs are always AF, similar to H-ZGNRs and 
ZGNRs without H-passivation V5 . The behavior of the AE as 
a function of the ribbon width is shown in Fig. [2] (a). For 
Nz > 4 the Aiip^ AF decreases almost linearly with increas- 
ing ribbon widths, in agreement with a recent calculation by 
Ong et alA However, in the case of 2-Fe-ZGNR the &E™$_ AF 
is much lower than that of 3-Fe-ZGNR. 

Contrary to the case of Fe-ZGNRs, we find that the in- 
teredge magnetic coupling of Co-ZGNR with a finite width 
is always FM, as shown in Fig. [2] (a). The absolute value 
of difference in total energy also decay as N z increases. It 
becomes negligible when the two edges are separated by a 
large ribbon width. By fitting the AE^ 1 ^ variation with w, 
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FIG. 2: (Color online). Total energy difference between the FM and 
AF states as a function of the ribbon width w. (a) Fe, Co and Ni 
terminated ZGNRs; (b) Fe and Co terminated AGNRs. Insets: width 
dependence of the exchange-correlation (XC) contribution as well as 
kinetic and electrostatic (KE) contribution to the total energy differ- 
ence between the FM and AF states. 



one obtains a decay law close to w~ 2 ' 5 . However, the almost 
monotonic width dependence of coupling for Fe-ZGNR and 
Co-ZGNRs shown in Fig.|2](a) are only part of the story and 
will be re-examined later. The smaller amplitude of interedge 
magnetic interaction in Co-ZGNR at large w as compared to 
Fe-ZGNR is primarily due to the fact that magnetic moments 
of Co-ZGNRs are more localized localized at the edges^. In- 
side the Co-ZGNR the magnetic tails of the edge moments 
decay much faster than Fe-ZGNR, resulting almost zero mo- 
ments in the inner C atoms of the ribbon. This localization 
effect is more remarkable in the case of Ni-ZGNRs. The in- 
teredge exchange interaction is found to be negligible in Ni- 
ZGNRs, as the FM and AF states become almost degenerate 
at even the shortest width with Nz = 2. 

The total energy difference between FM and AF states 
can be separated into non-interaction one-electron (kinetic), 
Hartree (electrostatic) and exchange-correlation (XC) contri- 
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A close inspection of the AE F ^_ AF curves in Fig|2](a) shows 
that the XC contribution to the total energy, is clearly respon- 
sible for the different type of interedge coupling in Fe-ZGNRs 
and Co-ZGNRs. For Fe-ZGNRs the AE™ _ AF is always pos- 
itive, while in Co-ZGNRs the AE F ^_ AF is negative but be- 
come negligible for large interedge distances. A compari- 
son between AE F ^_ AF and AE F ff_ AF indicates that FM cou- 
pling always lower the kinetic and electrostatic (KE) parts of 
the contribution for both Fe-ZGNRs and Co-ZGNRs, espe- 
cially for ribbons with short interedge distances. However, the 
two edges of Fe-ZGNR are still AF coupled resulting from a 
large energy difference from the XC contribution. Take the 
2-Fe-ZGNR as an example: the FM coupling lowers the KE 
contributions to the total energy by 41 meV/unit-cell but at 
a cost of -61 meV/unit-cell in the XC contribution, result- 
ing a AE F ^_ AF = 20 meV/unit-cell for an AF coupling as the 
ground state. 

For GNRs with armchair edges terminated by Fe and Co 
atoms, the interedge coupling between the two edges favor 
either AFM or FM depending on the ribbon width. The en- 
ergy difference between the FM and AF states as a function of 
the ribbon width is plotted in Fig|2](b). A damped oscillatory 
behavior, as shown in Fig|2](b) for both Fe-AGNRs and Co- 
AGNRs, clearly indicate a Ruderman-Kittel-Kasuya-Yosida 
(RKKY)-typeii interedge exchange. A 7V fl -Fe-AGNR is FM 
coupled only if N a - 3m + 2 (where m is a positive integer), 
and a 7V fl -Co-AGNR is FM only if N a = 3m. 

The well-known long-range RKKY interaction between 
two magnetic impurities in a non-magnetic host material is 
mediated by the conduction electrons of the host, and the cou- 
pling strength J can be written a a 15,16 



J(w) = J cos(2£ F w + (f>)/r D , 



(2) 



where D is the assumed dimensionality, w is the distance be- 
tween two impurities and kp is the wavevector at the Fermi 
level. By choosing dimensionality D = 1 and set w as the in- 
teredge distance of the ribbon, our results of the total energy 
difference AE F ff_ AF (w) °c J(w) can be fitted very well using 
Eq. (2) with a period = n/kp of about 3.7 A for both 
Fe-AGNRs and Co-AGNRs and slightly different phase c*. 
This is analogous to the case of the interlayer exchange cou- 
pling between ferromagnetic layers separated by nonmang- 
netic metallic spacer 17 , where a two-dimensional range func- 
tion, J(d) = J[)Cos(qpd + 4>)/d 2 , is used to describe the oscil- 
latory behavior of the coupling strength as a function of the 
spacer thickness d. The qp, which determines the oscillation 
period, is the critical spanning vector parallel to the interface 
normal that connects two sheets of the Fermi surface of the 
spacer at a point where they are parallel to each other. Simi- 
larly, here we find that the fitted 2kp = qp is exactly the criti- 
cal spanning vector connecting two inequivalent Dirac points 
K and K' in the Brillouin zone of graphene in the direction 
of the metal- AGNR interface (along the armchair edges) nor- 
mal, which is shown in Figj3] Another spanning vector q'p, 
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FIG. 3: (Color online), (a) Real space structure of the honeycomb 
lattice with two sublattices A and B in red (light) and black (dark) 
colors respectively, (b) The critical spanning vectors q F and q'js for 
graphene in reciprocal space. The q<p connects two inequivalent Dirac 
point K and K' and parallel to the k y axis which is the metal-AGNR 
interface (along x axis in real space) normal. The q'p connects two 
equivalent Dirac point K and K (or K' and K') and parallel to the k x 
axis which is the metal-ZGNR interface (along y axis in real space) 
normal. 

which connects two equivalent Dirac points K-K or K'-K', 
determines the oscillation period Tz of the interedge coupling 
for ZGNRs. However, following the discussion above we get 
Tz = 2n/q'p = 2.1 A which coincides with the interedge lat- 
tice spacing in ZGNRs. This clearly explains why the damped 
oscillatory behavior is not shown in actual width dependence 
of the A£^ al _ AF for either Fe-ZGNRs or Co-ZGNRs but a 
monotonic behavior. 

The width dependences of the XC and KE contributions to 



the total energy difference are also shown in Fig|2](b) for Fe- 
AGNRs and Co-AGNRs. Both AE™. .„ and AE%?. 4C show 

FM-AF FM-AF 

similar oscillatory behavior as the AE l °^_ AF with a same pe- 
riod, a phase shift as compared to the AE^_ Ap curve is found 
for the KE contribution in Fe-AGNRs, and for the XC contri- 
bution in Co-AGNRs. This indicates that the damped oscilla- 
tory behavior in the width dependence of AE l °^ AF is result- 
ing from a competition between different contributions to the 
total energy. 

In summary, we have presented a study of the interedge 
magnetic coupling of Fe, Co and Ni terminated graphene 
nanoribbons through first-principles calculations. We find that 
the ferromagnetic ordering of the metal terminations at each 
edge of the ribbon is favored for both ZGNRs and AGNRs. 
Whether the interedge magnetic coupling is ferromagnetic or 
antiferromagnetic depends to a large extent on the type of 
metal atoms and edges, as well as the ribbon width. The 
two edges for Fe-ZGNR are found to be antiferromagneti- 
cally coupled while for Co-ZGNR ferromagnetic coupling is 
favored, and the strength of the interedge exchange interac- 
tion decreases as the ribbon width increases. For both Fe- 
AGNRs and Co-AGNRs the interedge exchange interactions 
show damped oscillatory behavior as a function of the rib- 
bon width with a period of about 3.7 A. The interedge mag- 
netic coupling is negligible in Ni-ZGNRs, and Ni-AGNRs are 
found to be nonmagnetic. 
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